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Abstract 

FDTD  modelling  of  dielectric  PBG  structures  is  considered.  In  this  paper  PBGs  are  as- 
sumed to  be  related  to  layered  structures  which  are  especially  used  in  optical  waveguiding 
applications.  Many  essential  features  can  be  modelled  2-dimensionally. 

1.  Introduction 

Photonic  bandgap  (PBG)  material  is  periodically  inhomogeneous  material,  in  which  electromag- 
netic waves  having  frequency  in  the  bandgap  range,  can’t  propagate.  Due  to  the  periodical 
structure,  PBGs  are  often  called  photonic  crystals,  too  [1].  An  example  of  2-D  PBG  material  is  a 
uniform  lattice  of  cylindrical  air  holes  inside  dielectric  host  material.  In  real  3-D  world  this  might 
correspond  to  a silicon  plate  having  cylindrical  air  holes.  Properly  manufactured  PBG  material 
can  be  used  as  a frequency  selective  reflective  surface.  Potential  applications  are  highly  efficient 
optical  lasers  and  sharp  bends  in  optical  waveguides,  for  example.  In  this  paper  purely  dielectric 
PBG  structures  are  considered,  even  though  a PBG  can  also  be  manufactured  using  metal  rods. 
The  chosen  PBG  lattice  is  triangular,  formed  by  cylindrical  air  holes  in  dielectric  host  medium 
(figure  1).  The  main  interest  is  in  2-D  modelling  of  the  PBG  waveguides  and  waveguide  bends. 

The  behaviour  of  EM  fields  in  space-time  in  a PBG  structure  can  be  modelled  rising  FDTD 
(Finite-Difference  Time-Domain)  which  essentially  means  solving  Maxwell’s  equations  in  dis- 
cretized space  and  time  co-ordinates  [2].  FDTD  has  been  widely  used  in  PBG  research  during 
the  last  few  years,  see  [3],  [4],  for  example. 

2.  Modelling  PBG  Waveguides 

Modelling  PBG  structures  requires  some  important  properties  from  the  FDTD  software  used: 

• model  for  a dielectric  cylinder 

• freedom  to  describe  complicated  excitation  fields  i.e.  the  time-domain  behaviour  and  the 
field  distribution  in  space. 

• possibility  to  store  “measured”  fields,  at  certain  points  and  time  intervals,  to  hard  disk. 
Stored  time-dependent  fields  are  analyzed  later  (F-transforms,  power  flow  analysis  etc.) 

The  modelling  scheme  in  general  can  be  considered  as  a “virtual  measurement  setup”.  There 
must  be  a source  i.e.  a field  excitation,  for  which  time  and  spatial  dependence  is  properly  chosen. 
The  space- time  behaviour  of  EM  fields  is  governed  by  Maxwell’s  equations  which  are  solved  using 
FDTD.  Also,  there  must  be  measurement  surfaces,  on  which  time-dependent  fields  are  sampled 
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and  stored  to  hard  disk  for  later  analysis.  Figure  1 shows  measurement  setup  involving  straight 
PBG  waveguide. 

When  analyzing  the  performance  of  a guiding  structure,  computing  power  flow,  instead  of 
just  measuring  field  values,  should  be  preferred.  Just  inspecting  field  strenght  does  not  tell,  how 
well  the  power  is  actually  propagating.  For  example,  with  a standing  or  evanescent  wave,  field 
value  near  the  output  can  be  ^ 0 but  power  is  not  propagating. 

Power  flow  is  computed  by  spatially  integrating  the  normal  poynting  vector  S.  The  in- 
tegration is  done  over  a surface  (3-D  simulation)  or  over  a line  (2-D).  Using  simple  and  fast 
trapezoidal  rule  integration  leads  to  weighted  sum  of  poynting  vector  values.  For  example,  in 
2-D,  when  computing  power  flow  in  -hr-direction  through  measurement  line  Ml  (®  = x\,  y € 
[yi.lfe],  see  figure  1),  one  uses  stored  fields  Ev{xy,y,t)  and  Hz(xi,y,t).  Poynting  Sx(xi,y,t)  = 
Ey(xi,y,t)Hz(xuy,t)  and  time- dependent  power  flow  (per-unit-length)  is 


V2-yi  = {N-  l)Ay,  (1) 


where  the  approximation  corresponds  to  trapezoidal  rule  integration,  and  Ay  is  FDTD  cell-size. 

By  computing  power  flow  in  frequency  domain,  P(/),  one  can  effectively  obtain,  for  example, 
power  transmission  T(/)  = P^iD/Piif)  for  a certain  waveguide  bend.  This  requires  stored 
fields  Ei(r{,t)  and  i = 1...M,  where  M is  the  number  of  measurement  planes  (often 

M = 2).  Performing  Fourier  transform  from  time  to  frequency  domain  one  gets  Ei(ri}  f)  and 
Hi(ri,f)-  Then  for  every  point  r*  (essential  field  components  of  E and  H depend  on  the  plane  i 
orientation), 

S{ri,f)  — ^R[15(rj,/)  JT(n,f)]  (2) 

is  computed.  Via  integration  of  S,  one  gets  power  flow  through  measurement  planes,  Pj(/), 
i = Power  flow  analysis  in  frequency  domain  is  efficient,  because  only  one  simulation, 

using  pulse-like  excitation,  is  needed  to  get  the  transmission  T(f). 

The  absolute  values  are  not  needed  for  power  flow  at  different  frequencies  or  locations,  because 
only  the  spectrum  of  power  flow  matters  or  T(f)  = P2(/)/Pl(/),  for  example.  Hence,  formulas 
such  as  1 and  2,  are  used  without  coefficients  Ay  and  1/2.  Such  simplifications  are  used  also  with 
3-D  analysis.  Also,  the  P(f)-functions  are  divided  by  the  power  spectrum  of  the  input  signal,  to 
actually  see  the  filtering  effect  caused  by  the  waveguide  itself. 

In  the  formulas  above  it  was  quietly  assumed  that  E,  H and  5 sharing  the  same  co-ordinate 
(x,y,t)  are  also  physically  at  the  same  point  in  space-time,  which  is  not  exactly  true  due  to 
Yee’s  FDTD  cell  [2].  Additionally  to  numerical  dispersion  error,  also  this  may  cause  some  error 
at  higher  frequency  range. 


3.  Numerical  Examples 

The  numerical  examples  are  done  in  2-D  and  assuming  TEX  polarisation  i.e.  H = H(x,  y,  t)uz 
and  E = Ex(x,  y,  t)ux  + Ey(x,  y,  t)uy. 

First  example  is  a power  flow  analysis  of  a straight  PBG  waveguide  (linear  lattice  defect). 
Figure  2 shows  a ifs- field  snapshot  at  one  time  instant  and  figure  3 (left)  P(/)  at  one  measure- 
ment plane,  with  waveguide  width  w as  a varying  parameter.  PBG  material  remains  same  i.e. 
frequency  bandgap  is  same  in  every  case.  When  / belongs  to  bandgap,  PBG  works  as  a reflecting 
surface  or  as  a boundary  as  far  as  energy  flow  is  consented.  A mode  can’t  propagate  power  in 
the  waveguide,  if  / < fcut-off-  Cut-off  frequency  goes  up  as  the  waveguide  is  made  narrower. 
This  effect  can  be  seen  in  figure  3.  It  may  be  practical  solution  to  determine  the  rhannpl  width 
so  that  only  fundamental  mode  propagates  in  the  bandgap  frequency  region. 
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Second  example  deals  with  determining  the  wavelength  A (/)  of  a propagating  wave,  in  a 
straight  PBG  waveguide  which  can  be  approximately  considered  as  a lossless  periodic  waveg- 
uide, when  / is  in  the  bandgap  (small  radiation  loss).  As  in  example  1,  narrow-band  pulse-like 
excitation  was  used,  and  F-transformed  fields  at  two  planes  were  computed.  Using  the  field 
phase-difference  between  the  planes  (plane  separation-' was  a i.e.  the  waveguide  period),  wave- 
length A (/)  was  obtained.  Figure  3 (right)  shows  the  plot  A(/),  again  is  as  a parameter.  The 
results  were  checked  at  some  frequencies  using  time-harmonic  excitation. 

Note  that  in  examples  1 and  2 a PEC-boundary  at  y = 0 (figure  2)  was  used  to  exclude 
the  antisymmetric  modes.  So,  when  / is  in  the  bandgap  and  small  enough,  the  field  roughly 
corresponds  to  fundamental  mode.  Willing  to  investigate  antisymmetric  modes,  one  should  use 
PMC-boundary. 

Third  example  is  a transmission  analysis  of  a 60°  waveguide  bend.  P(f)  was  computed  at 
two  planes.  The  J3>field  snapshot  and  transmission  T(/)  are  shown  in  figure  4. 

Due  to  some  unwanted  reflections  the  field  sampling  time  window  must  often  be  limited.  This 
causes  some  decrease  in  accuracy  in  F-transformed  fields,  P(f)  and  T(f).  Considering  practical 
modelling  and  design,  a complicated  structure  is  much  more  easily  descibed  using  integer-like 
dimensions.  The  results,  a frequency  response  etc.,  can  be  easily  transformed  later  into  the 
frequency  range  of  interest  which  corresponds  to  the  real-world  physical  dimensions. 

4.  Conclusion 

2-D  modelling  of  PBG  structures  was  demonstrated  by  way  of  few  examples.  Using  pulse-like 
excitation  and  after  the  simulation  F-transforming  the  fields  at  measurement  planes,  one  can 
effectively  study  the  features  of  a PBG  structure.  It  seems  that  2-D  simulation  is  often  suitable 
for  design  purposes,  when  dealing  with  layered  structures  where  the  field  is  mostly  consentrated 
in  one  layer  (PBG  plate).  Thking  the  third  dimension  into  account  often  means  just  higW 
radiation  loss  and  a bit  frequency-shifted  response,  due  to  a different  effective  refraction  index 
seen  by  the  wave. 


Figure  1:  Triangular  lattice  PBG.  Straight  waveguide  measurement  setup. 


Figure  2:  If*-field  snapshot.  Computation 
domain  was  halved  using  PEC-boundary 
at  y — 0.  Consequently  only  symmetric 
modes  can  exist.  The  excitation  field  was 
chosen  to  support  1.  symm.  mode  i.e.  the 
fundamental  mode. 


Figure  4:  Left:  infield  snapshot.  Right:  power  transmission  of  a 60°  bend. 
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